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Abstract. - We generalize the concept of geometrical resonance to perturbed sine-Gordon, 
Nonlinear Schrodinger and Complex Ginzburg-Landau equations. Using this theory we can 
control different dynamical patterns. For instance, we can stabilize breathers and oscillatory 
patterns of large amplitudes successfully avoiding chaos. On the other hand, this method can 
be used to suppress spatiotemporal chaos and turbulence in systems where these phenomena 
are already present. This method can be generalized to even more general spatiotemporal 
systems. 



Spatiotemporal chaos PHEIHE] is one of the most important (and most studied) 
phenomena of recent years. Chaos can be advantageous in some situations, while in many 
other situations, it should be avoided or controlled [5,7 8 9 10 lllllll2lll3lll4| . In certain cases, 
the desired effect is a high-amplitude periodic oscillation. We should drive a nonlinear system 
with a large external force to produce such a high-amplitude oscillation. However, this should 
be done in such a way that chaos is avoided. Different feedback mechanisms have been devised 
to control chaos [7111511161117) . A great deal of research has been dedicated also to the problem 
of suppresing chaos by harmonic (or just periodic) perturbations |llH14l lT8. 19, 20 2"Tl l22ll2"31 
Among those works are the ones that use the concept of Geometrical Resonance 

(GR) [Il[ni[23El[23[13[27J. 

In Ref. 52 the concept of GR was used as a chaos-eliminating mechanism for the perturbed 
ip 4 equation. In this letter we generalize the concept of Geometrical Resonance to a very gen- 
eral class of spatiotemporal systems which includes the sine-Gordon, Nonlinear Schrodinger, 
Boussinesq, Toda lattice and Complex Ginzburg-Landau equations (among others). We will 
use this concept as a method of chaos control when these equations are nonintegrable because 
of the presence of perturbations. GR is an extension of the linear notion of resonance to a 
nonlinear formulation based on a local energy conservation requirement |23| . 
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Let us consider the partial differential equation 

K [4»}+K 1 [<f>,x]=q(x,t)P[<f>], (1) 

where Kq [<f>] and P [<fi] are functions of cf>, and its derivatives: <f>t, 4>x, 4>tt, <t>xx, etc. Equation 
Kq [(f)] — is an integrable Hamiltonian system. This can be, for instance, the sine-Gordon 
equation (SGE), the Nonlinear Schrodinger equation (NLSE), the Boussinesq equation, or the 
Toda lattice |28ll29j . On the other hand, K\ [cj>, x] includes dissipative terms and q (x, t) P [(f)] 
is a very general driving force 5 , 30 31, 32] - 

At GR the amplitude, frequency, and space-time shape of q (x, t) must satisfy some con- 
ditions so that some dynamical properties of the conservative system are preserved. We will 
call 4>Qn(x,t) a GR solution of Eq. if 

K 1 [cj>GR,x}=q(x,t)P[<f> GR }. (2) 

This implies a local energy conservation requirement. The energy integral that is conserved 
for equation Kq [(f)] — is locally conserved for Eq. Q if condition @ holds. We can use 
this condition as a mechanism for chaos control when an additional condition holds: the GR 
solution must be an asymptotically stable solution of the (full) equation Q). This condition 
is introduced here for the first time. We will call Eq. the exact GR condition and the 
solutions that satisfy this condition will be called GR solutions. 

We can consider the energy of the system as a "local almost adiabatic invariant" |33j . 
Then we can write an approximate GR condition 

dH\ 

where H is the energy of the system and T" is the period of the chosen solution of equation 
K [4>] =0. 

As an example, let us investigate the well-known driven and damped SGE 

4>tt - 4>xx + l4>t + sin (f> = q{x,t). (4) 

Suppose the task is to produce breathers of large amplitudes without entering a chaotic 
regime. The exact breather solution to the unperturbed SGE is 



<f) (x, t) = 4 arctan 



\/l — lo 2 sin (cut) 
o cosh (Vl — u) 2 x) 



(5) 



where lo is arbitrary in the interval lo 2 < 1. 

The external force q(x, t) satisfies the GR condition © when 



4 7 \/l - uj 2 cos (wt) 

q(x,t) = qGR (x,t) 



cosh (vT — U! 2 Q 



sin 2 (ujt) 
cosh(\/l—uj 2 x) 



(6) 



In Eq. Q if q(x,t) is given by I©, the function (0 is an exact solution of the complete 
Eq. (0J. When the parameters that define the perturbation I© are fixed, there is only 
one frequency for which function (|SJ| is the solution. This frequency is determined by that 
appearing in ©. 
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Fig. 1 - Robustness of the breather in Eq. (J1J. (a) Even with random initial conditions the breather 
is reorganized (u> = 0.707, 7 = 0.45). In the numerical simulations with the discretized equation the 
initial conditions were produced by a pseudorandom number generator of uniformly distributed values 
in the interval [—1,1]. (b) Irregular dynamics produced with a perturbation where the amplitude A, 
the frequency ui and the range Q are not close to satisfy the GR condition (A = 4.5, uj = 0.707, 
= 6.6). 

It is not difficult to show that the solution J5j is asymptotically stable in the framework 
of the full Eq. (@J with q(x,t) given by Eq. ©. In the framework of the unperturbed SGE 
breathers form a continuum of solutions similar to the periodic solutions around the fixed 
points called centers in Dynamical Systems theory. These solutions are stable in the sense of 
Lyapunov but they are not asymptotically stable. However, the breather solution JSJ in the 
framework of Eq. (0J with q(x, t) given by Eq. © is a spatiotemporal limit cycle. That is, 
this is a spatiotemporal attractor. All close initial conditions (in all space configurations) for 
t — > 00 will tend to behave as this solution. This phenomenon is shown in Fig. ^a). 

Perturbation qcn {x, t) can be approximated by a function of type q(x, t) = f(t)g(x) where 
g(x) is a bell-shaped function and f(t) is a time-periodic function. This kind of perturbations 
has been used in several studies of the SGE |M105] , 

The general study of Eq. using the GR concept and the breather solutions leads 
to the following conclusions: We can avoid chaos with amplitudes A of q(x, t) for which 
\A\ < 47VI — and lu 2 < 1. On the other hand, the range Q of the function g{x) (i.e. the 
interval of x where g(x) is not exponentially small) should be Q < -^==^. 

In some cases, when these conditions are not satisfied, the breather is not stabilized and 
we can get an irregular behavior. Figure^b) shows an example of the dynamics produced by 

Eq. Q with q {x, t) = A cos (ut) cosh | + (^) wher e A = 4.5, u = 0.707, 

Q = 6.6. 

There is a wealth of works |34ll55| dedicated to the numerical investigation of perturbed 
sine-Gordon equations using external forces of type q(x,t) = f(t)g(x). All the results are in 
agreement with our theoretical results. 

Sometimes we have the task of suppressing chaos using a given harmonic perturbation. 
Consider e.g. the following equation: 

0tt - 4>xx + 70t + sin = /o sin (uj d t) + f c sin [u c t + 9} , (7) 



where / sin (w^t) is a chaos-producing excitation, while f c sin \o c t + 9] is a chaos-suppressing 
excitation. The parameters of the chaos-suppressing excitation should be determined. In this 
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Fig. 2 - Suppression of spatiotemporal chaos using a given chaos-suppressing excitation in Eq. @ 
(/o = 0.91, u>d = 0.6). (a) Well-developed spatiotemporal chaos when f c — 0. (b) Controlled dynamics 
with f c = 0.4, uj c = 0.6, 9 = tt/2. Controlling spatiotemporal chaos using a localized excitation in 
Eq. (c) Spatial profile for a given time moment corresponding to spatiotemporal chaos when 

F c (x,t) — 0. (d)-(p) Spatial profiles for different time instants when F c (x,t) is given by Eq. JB), 
7 = 0.1, fo = 0.5, Lo d = 0.6, w = 0.6. The time instants are (d) t = 0.97, (e) t = 7.9, (f) t = 15.4, (g) 
t = 21, (h) t = 26.2, (i) t = 36.8, (j) t = 42, (k) t = 46.8. 



case we can use the condition (0 to find parameters for the chaos-suppressing perturbation. 
Fig. IU(a)-(b) shows an example of chaos control using this technique. 

We would like to stress here that the force © can be used to control a well-developed 
spatiotemporal chaos. 

Let us consider the following equation: 

(f>tt - <t>xx + l4>t + sin (j) = f sin (u d t) +F c (x,t). (8) 

When F c (x, t) = 0, the system presents spatiotemporal chaos for — oo < x < oo, see Fig. 
I2jc). Now, if we turn on the controlling force F c (x, t) defined as function JBJ, we obtain a very 
regular spatiotemporal pattern as that shown in Fig. |3(d)-(k). The most important remark 
here is that we are controlling spatiotemporal chaos in the whole space using a localized 
perturbation. We should add here that other works have used localized perturbations to 
control spatiotemporal chaos |36U37ll55] , 

The damped and ac-driven NLSE 

icbt + cbxx + Zlcbfcb + iacjj^ee^ (9) 

is another fundamental model in many areas of physics |39ll40ll4"T] . At sufficiently large e the 
dynamics of this model becomes chaotic |41| . 

Suppose we have a general driving term: i<fi t + <fi X x + 2 |0| 2 <f> + = Q ( x , t) ■ We will 
take the one-soliton solution of unperturbed NLSE [2H1 as a GR solution: Then q(x,t) must 
satisfy the condition: 



q GR (x, t) = — 

COSh (y/OJX) 



(10) 



J. A. Gonzalez, A. Bellom'n, L. I. Reyes, C. Vasquez and L. E. Guerrero: Geometrical Resonance in Spatiotemporal Sy 



where u> can be any positive number. Thus, if the perturbation is localized and the amplitude 
e satisfies the condition e w a^/uj, then the chaotic regime can be avoided. The one-soliton 
solution of NLSE is stabilized. We should remark that this can be achieved also by other 
localized perturbations. 

If we use the two-soliton breather solution as a GR solution, we can obtain another driving 
force satisfying a GR condition: 

4a [cosh (3x) + 3e j8t cosh (x)] e^ t+ ^ 
q G R[x,t)- cosh ( 4a ,) + 4cosn ( 2x ) + 3 cos (8t) ' (H) 

See Ref. [29 for a discussion of multisoliton solutions. As in Eq. I|1U|) we can introduce here 
an arbitrary parameter u>. However, in this case, we are more interested in the relationship 
between the two intrinsic frequencies of the solution. This force can be approximated by a 
function of type 

F(x,t) = e 1 g 1 {x)e i (» lt+ "M + e 2 g 2 {x)e l{ ^ t+v ' 2 \ (12) 

where u>i = 1, u>2 = 9, and <?i(x) and g 2 {x) are localized functions. 

In Ref. a breather was stabilized using a two-frequency drive: F (x,t) = E\e %u}lt + 
e 2 e 1 " 2 *, where w\ = 1 and lo 2 = 9. This result can be seen as a confirmation of the GR 
approach for the NLSE. We should add that the phenomenon of breather stabilization is quite 
robust. For instance, if u>\ = 1, in addition to oj 2 = 9, other close frequencies can be used, 
namely lo 2 — 8 and oj 2 — 10. This means that the GR condition can be satisfied approximately 
and that the Eq. (|3J) can also be used as a guide for the search of a controlling force. As 
in the case of the breather of the SGE, the breather of the NLSE is asymptotically stable. 
We have checked numerically that the nonchaotic breather solution is the most stable one. 
Unfortunately we have no space to show a picture. 

Regarding localized excitations we should emphasize that GR analysis explains diverse 
fundamental results on stability of localized solutions previously obtained by perturbation 
theory |14ll5ll41j including those relative to one-soliton and two-soliton solutions of the NLSE. 
In this sense, in future works, it would be interesting to consider the case of the N-soliton 
solutions (see e.g. Appendix B, Ref. [2*^)1. 

The control of spatiotemporal chaos (or turbulence) in the Complex Ginzburg-Landau 
equation (CGLE) 12 , 42 : 43 : 44 : 45 : 46 , 47 is a problem of great practical interest We are 
interested in the modified CGLE [Hll^ : 

</H=4> + (l + ic x ) (ft xx - (1 - «c 3 ) |0| 2 (ft + F c (x, t) . (13) 

The term F c (x, t) is the control signal. Without the control signal (F c (x,t) = 0), the 
turbulence develops when the Benjamin- Feir condition 1 — C1C3 < is satisfied. This equation 
can be rewritten in the following form: 

i(ftt + cicftxx + C3 \<ft\ 2 (ft = i \ 4>xx +<ft- \<ft\ 2 <ftj + iF c (x, t) . (14) 

When the right hand side of Eq. (|14|) is zero, it reduces to the NLSE. 

If (ft (x, t) = f (x) exp (icot) is a soliton solution of the NLSE, then we can use the following 
controlling signal 

F c (x,t)=[f(x)-f(x)-f xx (x)]exp(iut). (15) 

Equation (|13[) (with F c = 0) presents turbulence for c\ = 2, C3 = 0.8. We have been able 
to suppress this turbulence using the F c (x, t) given by Eq. 1)15(1 with to = 12 and / (x) is the 
one-soliton solution of equation cif xx — ujf + c 3 / 3 = [2"§| . 
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In this context, we should explain that in some cases, the stabilization process can require 
a force that is not a small perturbation. Furthermore, this technique can be used both as a 
way to stabilize a pre-existing solution of the unperturbed system and as a way to impose 
an arbitrary solution to the system. However, the success of all these endeavors depends on 
a very important fact: the final solution should be an asymptotical stable solution of the 
perturbed system. Incidentally, we should mention that the stabilization of unstable plane 
waves in the CGLE can be done using a nonlinear diffusion term |48|. 

In some situations we can apply some perturbations using technological means in order 
to satisfy the stability conditions. Nevertheless, we should say that, very often, nature itself 
can apply the controlling perturbations. There are many natural regimes described by the 
mentioned equations in the presence of perturbations where the resulting dynamics is not 
chaotic. Our results can provide an explanation for these phenomena. 

Numerous observations and experiments show that elastic waves from natural phenomena 
and human-made machines may alter water and oil production |49j . In some cases wave 
excitation may appreciably increase the mobility of these fluids. A new technology |49j based 
on these experiments is used to stimulate the reservoir as a whole. Here seismic frequency 
waves are applied at the earth's surface by arrays of vibrators. Many of the phenomena 
involved in this effect are described by the equations discussed in this paper, namely: NLSE, 
SGE, Boussinesq equation and other equations of type |QJ (see Refs. jSU]). For the optimization 
of the method, it is necessary to sustain spatiotemporal nonlinear oscillations of the reservoir 
with some frequency and shape. Based on ideas related to the results presented in this paper 
we have designed a new technology using a specific geometrical arrangement of the surface 
vibrators |51| . 

The nonlinear PDE possess an infinite number of different solutions. Among them one 
can choose a feasible one in order to implement. Even if only a given type of perturbation 
is allowed due to technical limitations, it is always possible to use the approximate condition 
© as in the case of task Q . 

The concept that links all the situations where we have been able to suppress chaos is 
based on the mutual cancellation of nonintegrable terms as described by equations J5J) and 
(JsJJ. In other words, we should add some temporal perturbation in such a way that (at 
least approximately) both the dissipative and the total driving terms mutually cancel. A 
remarkable situation (which is a particular case of the general theory but, at the same time, is 
present in all the studied systems) is that of breather-like oscillations. These patterns can be 
stabilized using some spatially localized time-periodic perturbations, where the amplitude, the 
spatial range and the frequency must satisfy some relationship. However, this phenomenon is 
robust. A fine-tuning is not necessary. There is always a whole valid interval of values for the 
amplitude, range and frequency that produces qualitatively the same result. 

The most common perturbation in scientific research is F (t) — /q cos (u>t). However, 
nature is very rich in dynamical behaviors. Our work shows that using very general spa- 
tiotemporal perturbations F (x, t) we can make the difference between regular or chaotic 
behavior. Using certain spatiotemporal perturbations F (x, t) we can stabilize a breather or 
we can produce a turbulent dynamics. We have been able to control different patterns in 
the sine-Gordon, Nonlinear Schrodinger, and Complex Ginzburg-Landau equations. Each of 
these systems possesses wide applications in many areas of Physics. Furthermore we believe 
that these ideas can be applied to other systems. 
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